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Performance Flow Reactors 


Various Levels 


CHOLETTE?, BLANCHET? and LEONCE 


Chemical flow reactors studied include models 
which partial mixing and piston flow are present and 
others which involve partial mixing and short-circuit. 
Mixing levels considered extend all the way between 
the extreme cases tubular and continuous-stirred 
tank reactors. Various orders reaction are taken 
into account. 


The curves presented, obtained through analytical 
methods, show either the variation conversion 
with mixing level given values residence time, 
relative residence times for given conversion. 


Graphical methods solution are proposed. They 
involve, part, means computation whereby 
the overall conversion for number reactors 
series obtained from the individual conversions. 


While isothermal conditions only have been con- 
sidered, more interesting results are expected when 
other conditions are applicable, whether isothermal 
not. 


procedures pertaining chemical flow reactors are 
readily available for relatively simple models such 
tubular continuous-stirred-tank reactors (CSTR). 

Flow patterns such reactors correspond limiting cases. 
While CSTR’s can easily realized, tubular reactors seldom 
meet perfectly the assumption true plug flow. 

Actual reactors will often perform such way that results 
lie between those predicted for tubular and CST 

Various methods approach have been suggested deal 
with departures from ideal conditions. Residence time studies 
have been the object considerable work but 
not account for complete solution. Different 
have been used such work, some involving radioactive 
Thus, some have compared actual reactor the equivalent 
number CSTR series giving the same results. 

More recent studies were aimed analyzing reactor per- 
formance through the introduction diffusivity, involving 
Peclet number 19, 20, 

efficiency mixing continuous flow systems, different 


received September 25; accepted November 20, 1959. 
2Université Laval, Quebec, Que. 
Contribution from the Chemical Engineering Department, Université 


Laval, Quebec, Que. and presented the Chemical Engineering 
Conference, Hamilton, Ont., November 9-11, 1959. 


Les auteurs étudient des réacteurs chimiques 
régime continu comprenant des trouve 
tantét partielle frontal, 
tantét partielle court-circuit. Ils con- 
sidérent des niveaux d’agitation variant entre les cas 
extrémes réacteurs tubulaires réacteurs 
parfaitement agités. Ils présentent les résultats 
rapportant différents ordres réaction. 

obtient par méthodes analytiques des courbes 
qui montrent, soit variation conversion avec 
niveau d’agitation pour temps séjour donné, 
soit les temps séjour relatifs pour une conversion 
donnée. 

Les auteurs proposent des méthodes solutions 
graphiques. Elles comportent, partie, une facon 
conversion totale, pour des réacteurs 
série, partir des conversions individuelles. 

considére dans présent travail que des 
conditions isothermiques, mais prévoit aussi des 
résultats intéressants pour d’autres conditions, iso- 
thermiques non. 


models have been presented which various combinations 
piston flow, effective volume mixing and short-circuit were 
discussed. The case chemical reactors was introduced and, 
indicated, was expanded. 

The present study deals with models which partial mixing 
and piston flow exist and others where partial mixing and short- 
circuiting take place. 

Before investigating such reactors, useful review 
briefly the extreme cases: the tubular reactor where mixing 
present and the CSTR which exhibits perfect mixing. 

Design procedures rely deeply upon reaction rates. These 
involve number variables and can expressed 
different ways One the main factors temperature, 
heat effects can play important role the design reactor. 

While only isothermal conditions have been considered 
the present study, the scope the latter not limited that 
case alone. Endothermic exothermic reactions resulting 
non-isothermal conditions can studied well and will 
the subject future publication. 


TUBULAR REACTORS 


Design equations available for tubular reactors, applied 
mode! such shown Figure give the following expression, 
when the density liquid system, the number moles 
gaseous system, assumed constant: 
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One Reactant 


the many ways which express the reaction rate, 
the usual one for homogeneous irreversible reactions 


where the order the reaction. the velocity rate 


For first-order reaction, integration (1) gives: 


Smith (MacMullin and Weber and Eldridge and 
used the symbol D), the feed and effluent concentrations 
any flow system are related follows: 


Cy Cz = (4) 


plot which appears Figure From this graph, the 
residence time conversion can obtained directly. 
For second-order reaction, and Equation (1) 


yields: 


When plotted Figure convenient for determining 
either the residence time the conversion. combinations 
are also possible. 


For third-order reactions, and the design equation 


(7) 
x,)? 


and can useful when evaluating the ratios 
6/6, V/V, dealing with partial mixing, has been done 
later sections. Integration Equation (1) gives the following 
general expression for all values except unity (in which 

case Equation (5) 


Two Reactants 


Second-order reactions are often those which result from 
the interaction two components. For reactions the type 


products 
the reaction rate 


For the study these reactions often convenient define 
the ratio 


Ay 


Ca, 


this case, the degree conversion refers the limiting 
reactant, that lesser concentration. 


When the case reverts one second-order 
reaction with one reactant, and Equation (6) applies. 

Curves showing the variation C4, with different 
values are presented Caddell and Hurt and 
Lessels 

These authors, well Jenney and Nord also 
present nomographs relating Equations (5), (6), (7) and (9). 
While they serve similar purpose Figures they 
not apply the same way and are less convenient for some 
phases the present work. 

When deriving the preceding expressions for tubular 
reactors, true plug flow assumed. practice, such conditions 
are seldom attained. study the effect turbulent and 
laminar flow reactors was presented Bosworth 
Denbigh reports that while turbulent flow approximates 
closely the hypothesis true plug flow, laminar flow does not. 

indicated the latter and also Greenhalgh 
and Gorrigan and Young obtained for tubular 
reactors are also applicable batch reactors when the residence 
time the former the same the reaction time the latter. 

tubular reactor can considered made infinite 
number infinitely batch reactors continuous-flow 
stirred-tank reactors series the same total volume. 


CONTINUOUS-FLOW STIRRED-TANK REACTORS 
(CSTR) 


composition and temperature throughout the reactor. The 
concentration the outlet stream then the same that 
the reactor. Other pertinent assumptions are applicable. 

Many design equations are available, but the basic CSTR 
equation may written follows 


One Reactant 

general, the reaction rate for homogeneous, irreversible 
reaction given by: 

Some investigators have found convenient use 
this simple expression determine experimentally the velocity 
rate constant, 

Combining (10) and (11), one 

(a) for first-order reaction 


M 


(b) for second-order reaction 


(13) 
(c) for third order reaction 


While not presented here, graphs obtained from these three 
equations are similar form those presented for tubular 


Two Reactants 


the case second-order reaction the rate 
expressed follows: 
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The 


Integration Equation (1) gives this case: 
5 
C; 
2 
Equation (3) can thus written: 
: Fig 


Figure 


10° 
Figure determination variables tubular reactors, 
10° 10° 10" 10° 
Figure 3. Graphical determination of variables in tubular reactors, 
le second-order reactions. 


Figure 4 Graphical determination of variables in tubular reactors, 


14) j Figure 5— Continuous-stirred-tank reactor (CSTR). third-order reactions. 
lar are presented Lessels and Jenney 


‘The design equation obtained from (10) then; 


pointed out Reman tubular and CST reactors can 
equimolar concentrations, this equation reverts compared two ways: 
Equation (13). (a) the conversions obtained for given residence 
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10° 

> 4 

10" 


(b) the residence times volumes required achieve 
given conversion. 

The latter case has been the subject 

ies 11, 12, 26, 37) which not need extended for the 


present work. 
for the conversions obtained for given residence time, 


simple relation can derived from Equations (8), (10) and 
(11), the case one reactant, noting that 


1 x. 


which can rewritten the following form: 


1 
Xu 


This expression general one for all values except 
unity. this case, 


(17) 


Plots Equations (17) and (18) are shown Figure 
Some these have already been given Reman Similar 
curves can presented for any value those for 0.5, 
0.75 and 1.5 have been calculated but were not included 
Figure prevent overcrow ding. 

Figure for isothermal conditions. limited graph 
the inverse ratio has been presented Corrigan and 

While the ratio x,/xy close unity low and high 
values conversions, reaches maximum intermediate 
values. Although for the reactions under consideration the 
present work tubular reactor always more effective than 
CSTR equal residence time volume, more 
the region where the values x,/x through maximum 
the curves Figure may noticed that the particular 
value the conversion which the maximum occurs shifts 
with the order reaction. For each order reaction the 
maximum can read directly the curves Figure 
can calculated dividing Equation (17) xy, and setting 


d Xy 
Figure along with the corresponding values and xy. 

Since the conditions for which Figure applies are met only 
particular cases, another comparison the relative perfor- 
mance the two reactors can made showing the variaticn 
x,/xy function for given values the conversion 

may seen from this figure that the relative advantage 
the tubular reactor over the CSTR more pronounced 
intermediate conversions for the usual reaction orders, and 
Except low conversions, the tubular reactor always has 
decided advantage over the CSTR for first-order reactions. 


CONTINUOUS FLOW REACTORS WITH PARTIAL 
MIXING 


the models defined Cholette and Cloutier only 
the following are considered the present work: 

(A) Partial Mixing and Piston Flow. 

such model, fraction the reactor contents 
assumed perfectly mixed while the other, 
travels piston flow. Within this model, many combinations 
are possible depending the location the well-mixed zone. 
The limiting cases considered are described below and both 
analy tical and graphical methods are for their solu- 
tion. comparison the two models also made. 


(a) Model TM. 


shown Figure 10(a), the flow this model such 
that the reactants first undergo piston flow and then pass through 
zone perfect mixing. 


(b) Model MT. 


Figure 11(a) shows that the reactants enter zone perfect 


mixing and then pass into one piston flow. 

(B) Partial Mixing and 

The flow pattern and mixing characteristics are shown 
Figure 12. This model has been discussed briefly for zero- 
order reactions and for first-order reactions when short- 
circuiting present. studied more extensively and 
general way the present article. 


Partial Mixing and Piston Flow. 
Model 
Conversions obtained for given residence time. 


The flow through the reactor, Figure 10(a) such that the 
reactants first undergo tubular plug flow and then enter zone 


q 


perfect mixing, volume mV. The residence time, 


for the tubular section, assumed the same for all elements. 

Such model would behave like tubular reactor volume 


shown Figure 10(b). 
recent paper Zwietering 


Referring Figure and applying Equation (1), the results 


for the tubular part the reactor are: 


For the zone perfect mixing, applying Equation (10), gives: 


performance such reactor can then obtained 
combining Equations (19 and 20) and comparing the results 
those that would obtained with tubular reactor equal 
residence time, given Equation (1), being the theoretical 
outlet concentration for tubular flow. 
Combining Equations (1) and (20) gives: 


which, when introduced into Equation (19), leads the general 
expression: 


For first-order, irreversible, homogeneous reaction, Equation 
(22) simplifies to: 


Cr/Cr -( + m | 


or, terms conversion: 


Curves showing the variation with are given Figure 
13, values (V/q), corresponding anumber conversions 
the residence time increases, giving larger values 
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Figure 


Figure 


the 
Equati 
n= 


(Cr 


The 


! 
0. 
Xn 
4 a 
Ch 
0. 
( 
q 
(23) 


on 


TUBULAR 

t 


Figure 6—Conversion tubular reactor compared that CSTR 
equal residence time, for various orders reaction. 


1.40 


Figure 8- Variation of the maximum value of the ratio x,/x as a function 
the order reaction. 


Figure 10—Description Model TM. 


the conversion less and less affected the degree mixing, 
especially low values 

second-order reaction, when one reactant present, 
Equation (22) gives: 


m= 


1.40 


1.35 


1.25 


1.20 


1.05 


1.00 
0.0 0.2 0.4 0.6 0.8 1.0 


Figure 7—Variation of the relative conversion x; x,;, for reactors of 
equal residence time, as a function of the conversion in a tubular reactor. 


1.40 


Figure 9—Variation of the relative conversion x, x, as a function of the 
order of reaction, at given values of x. 


or, terms conversion: 


(26) 


Figure shows the variation with this case for 
given values 
two reactants unequivalent molar amount undergo 
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0.7. 

0.6 

a 
0.2 
0.1 

ne 

max 

1.30 0.4 

xX 

1.20 
0.95 
q 


Figure 11—Description Model MT. 


Figure 13—Variation the conversion with for first-order reactions, 
applicable Models and MT. Each curve corresponds residence 
time such that the conversion in a tubular reactor is read at m = O. 


Figure 14—. Variation of the conversion with m, for second-order reactions. 
Each set of curves corresponds to a residence time such that the 
conversion in a tubular reactor is read at m = O. The upper curve of 
each set for Model reactors, while the lower curve for Model MT. 


a=2 


MODEL 


Figure 16—Variation relative residence times with Model TM, 
for second-order reactions, at the conversions indicated on each curve. 


second-order reaction, the following expression derived from 
Equation (22 


(1 x) [x 1) + (R —_ x) In (aia) ] = 


(27) 
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STAGNANT ZONE 


Figure 12—Description of a reactor where partial mixing and short- 
circuit take place. 


MODEL 


0.70 


0.0 


0.74 0.78 0.82 


Figure 15—Variation with Model TM, for different orders 
reaction. The residence time such that the conversion tubular 
reactor would be 90%. 


MODEL 


Figure 17—Variation relative residence times with Model TM, 
for different orders reaction, fixed conversion 90%. 


might expected, examination this expression shows 
confirmed inspection Equations (9) and (15) for Oy. 

When the two reactants are present equimolar concentra- 
tions, and Equation (26) applies. 

graph can made from Equation (27), presenting families 
curves showing the variation with given values 
x,and Applying Equation (22) reactions order 0.5, 1.5 
and the following expressions are derived: 


sion, 

that 
conve 

sion 
dition 
less 
wouk 


mixin 


and 
equat 


under 


For 
t 
1.0 
> \ | | 
(10), 
0.9 
03 
1.0 
0.8 
0.3 1.0 


For 0.5 


For 1.5: 


For values and corresponding given conver- 
sion, obtained successive approximations. 

Application Equation (22) zero-order reaction shows 
that irrespective the values and the degree 
conversion depends only the residence time. 

For residence time such that the conversion tubular 
reactor would 0.9, figure (15) shows the variation conver- 
sion with for the above cases. may observed that con- 
ditions depart from those tubular reactor, the conversion 
less and less affected the extent mixing the order 
reaction increases from value 0.5; although curves are 
presented for very low orders reaction, the opposite effect 
would noticed varies from 0.5. 


Residence times volumes required for given 


For the tubular section the reactor shown Figure 
(10), the following equation holds: 


which, upon integration for gives: 


the reaction takes place thereafter zone perfect 
mixing, one obtains: 
C; Cr 
(33) 
which can transformed to: 


(34) 
Combining Equations (32) and (34) yields: 
valid for reaction orders other than one. 
plug flow prevailed throughout the above reactor, the feed 
and effluent concentrations would related the following 


The ratio residence times that the volumes 
obtained follows for identical conversions the two reactors 
under consideration: 


0, Eq. 36 
zero-order reaction, consideration Equations (35), 
(36) and (37) shows that for all values 
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For second-order reaction, Equations (35), (36) and (37) 
give: 


wherein different from zero unity. 

Curves showing the variation with different 
conversions, are given Figure the conversion 
increases, the relative residence times increase more and more 
rapidly when the degree mixing approaches unity. 

When two reactants undergoing second-order chemical 
reaction are present unequivalent molar concentrations, the 
following expression applies: 


The ratio then obtained follows: 


Many families curves can obtained from Equation (40), 
depending upon the values the different parameters contained. 
Only the particular case equimolar concentrations considered 
here, Equation (40) reverts Equation (38) when 
Curves for this particular case are shown Figure (16). 

When 0.5, Equations (35), (36) and (37) yield: 


. (413 
which must different from unity. 
For 1.5, (35) becomes: 


for Equation (35) gives: 


| 


Finding values Equations (42) and (43) approxima- 
tion and dividing the corresponding value Equation 
(36), gives the ratios for given conversion different 
values 


Such curves, well others for different orders reaction, 
are presented Figure (17), when the conversion 90%. 
Relative residence times are seen vary rather rapidly with 
for most the orders reaction, the degree mixing 
approaches unity. 


a 


oouw 


MODELS TM&MT 


Figure 18—-Variation of relative residence times with m, in Models TM 
and MT, for first-order reactions, at the conversion indicated on each 
curve. 


Figure 20—(a) reactor with a zone of perfect mixing located between two 
sections of tubular flow. (b) reactors in series, with a CSTR located 
between tubular reactors. (c) reactor with zone perfect mixing 
extending throughout the length of the reactor and surrounded by a 
section in tubular flow. (d) CSTR and tubular reactor in parallel. 


0.70 0.74 0.78 


0.82 


Figure 23 Comparison of the variation of conversion with m, in Models 
TM and MT, fora = 0.5, ata residence time such that x, = 0.9. 


For the particular case first-order reaction, Equations 
(31), (33) and (34) show that the feed and effluent concentrations 
the tubular section are related by: 


Cr j 
and those the well-mixed section 


MODEL 


0.70 0.74 0.78 0.90 


Figure 19—Variation of x with m, in Model MT, for different orders of 
reaction. The residence time is such that the conversion in a tubular 
reactor would 90%. 


a=2 
MODEL 


Figure 21—Variation of relative residence times with m, in Model MT, 
for second-order reactions, at the conversions indicated on each curve. 


1.0 


MODEL 
0.0 


Figure 22—Variation of relative residence times with m, in Model MT, 
for different orders reaction, given conversion 90%. 


0.0 


0.70 0.78 0.62 0.90 


Figure 24 -Comparison of the variation of conversion with m, in Models 
I'M and MT, for first and second-order reactions, at a residence time such 
that x, = 0.9, 
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Combining these two equations results 


Equation (5), seen that, for identical values 
conversion, 


The ratio for first-order reaction, thus 
obtained 


0.70 0.74 0.78 0.82 0.90 


Figure 25- Comparison of the variation of conversion with m, in Models 
TM and MT, fora = 1.5, at a residence time such that x, = 0.9. 


27—Comparison the variation relative residence times with 
m, in Models TM and MT, for a r 0.5 anda = 1, at a given conversion 
of 90%. 


Figure 29--Comparison of the variation of relative residence times with 
m, in Models TM and MT, for second-order reactions, at a given 
conversion of 90%. 
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the conversion approaches unity, the value increases 
rapidly and the term Equation (47) becomes 


The ratio residence times volumes for complete conver- 


sion simplifies to 
= 
= 1 


Curves for first-order reactions showing the variations 
6/0, with different values conversion are shown Figure 
18. The effect mixing relative residence times similar 
that shown Figure (16) for second-order reactions, except 
that more pronounced for first-order reactions. 


0.70 0.74 0.78 0.82 


Figure 26—Comparison of the variation of conversion with m, in Models 
TM and MT, fora = 3, at a residence time such that x, = 0.9. 


Figure 28—Comparison the variation relative residence times with 
m, in Models TM and MT, fora = 1.5, at a given conversion of 90°. 


Figure 30—Comparison of the variation of relative residence times with 
m, in Models TM and MT, for third-order reactions, at a given conversion 
of 90%. 
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Some the above expressions become indeterminate 
and the ratio more readily computed for all these 
cases means the following expression: 


Some expressions are also indeterminate can 


seen readily that the ratio 6/6, always equal unity this 
case because the reactor considered tubular one. 


Model 
Conversions obtained for given residence time. 


This model shows conditions similar those model 
except that here the reactants enter zone perfect mixing 
volume and then travel piston flow through section 
Figure Such combination has been discussed 
Zwietering and also Trambouze and Piret recent 
study dealing with the combination reactors giving optimum 
conversion. 

Referring Figure (11) and applying Equation (10) the 
feed and intermediate concentrations the reactors are related 
follows: 


For the tubular zone, Equation (1) gives: 


comparison similar that developed for model the 
reactor was tubular the inlet and effluent concentrations would 
(8). The latter, combined with Equation (52), for reaction 
orders other than unity, gives: 


Combining Equations (8) and (53), one obtains: 


which, when introduced Equation (54), gives the desired 
relations. 


For zero-order reactions, should expected, Equations 
(54) and (55) give the same results those obtained with model 


For second-order reaction, involving only one reactant, 
combination Equations (54) and (55) gives: 


Cr/C, 


or, terms conversion: 


Figure (14), values V/q corresponding various conver- 
sions, 

When two reactants are present, the following expression 


obtained: 


When the case reverts one second-order 
reaction and (57) applies. 

For other orders reaction, the combination Equations 
(54) and (55) gives the following expressions: 


fora 0.5 
= = J \ 
for 1.5 
fora 
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Values Equations (60) and (61) are obtained 

Curves showing the variation with for these cases are 
presented Figure (19) for residence time such that would 
90%. The same remarks apply this case for those 
Figure (15) pertaining Model TM. 

For the particular case first-order reactions, the application 
general Equations (52), (53) and (1) Model leads 
the same expressions, (23) and (24), were obtained 
for Model TM. The curves Figure (13) can then used 
this case. 

Indeed, for first-order reactions, can demonstrated that 
the results will always the same for given reactor, matter 
where the zone perfect mixing located, long the resi- 
dence time for all elements going through that zone the same, 
and equal Figures 20(a) and (b). would not apply 
reactor such shown Figures 20(c) and (d) because the 
residence time for elements going through the zone perfect 
mixing this case 


Residence times volumes required for given 
conversion. 


For the CSTR section the reactor, referring Figure 11, 
the following equation applies: 


the tubular section that follows, one obtains: 


For each reaction order, values are found for different 
given conversion, combining Equations (62) and (63). 
They are then divided the value obtained from 
(1) for the corresponding value The ratio thus 
obtained, function for given values the conversion. 

Here again, zero-order reaction gives for all 
values 

first-order reaction, since Model behaves like 
Model TM, Equations (4+) (50) inclusive are applicable, 
Figure 18. 

For second-order reaction, when the procedure 
described above leads the ratio: 


The ratio this case would thus be: 


6, 1 R- 

R-1 


Again, when the case reverts one second-order 
reaction and Equation (64) applies. 

For reaction order 0.5, when the ratio 
found be: 


When 1.5, the ratio involves approximation 
calculations. The first obtained from the expression: 


The value then calculated for the same conversion, 
using Equation (1), and the ratio easily obtained this 
manner. 


this case, the 


8(1 


(69) 


For conversion 90%, Figure shows the variation 
with for the different orders reaction discussed above. 
Unlike the performance model TM, seen that model 
the relative residence times vary more and more rapidly 
with the degree mixing the order reaction increases. 


2(1 


Curves derived from this expression are presented Figure 21. 
interesting note that for the limiting case complete 
conversion, when the ratio given Equation (64) 
reduces the expression 


(0/ 6.) = 


which happens the same that Equation (50) for first- 
order reactions. 


When two reactants are present and the following 
obtained: 


R — 


For purposes comparison, conversions have been calculat- 
for Models and residence time such that 
90%. The results, for different orders reaction, are presented 
Figures 23, 24, and and correspond the solid curves. 
may happen that given level mixing corresponding 
certain value the well-mixed section the reactor 
located between two zones tubular flow, shown Figures 
20(a) and (b). Depending the location this zone, the result- 
ing conversion will the shaded area horizontal line 
connecting the and curves. This valid provided the 
combined residence time between for all elements 
going through the tubular section. 

Inspection the preceding figures reveals that residence 
time studies, such suggested Sherwood ‘”, are not sufficient 
predict the conversion given reactor. necessary indeed 
know also the history the various elements they flow 
through the reactor since, given residence time, different 
conversions can obtained depending the pattern mixing. 


2 . 
23, 24, and given values For reaction orders 
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superior unity, Model gives higher conversion than 
Model MT, while the reverse true for values inferior 
unity. the case first-order reactions, equal conversions are 
obtained for both models. 

different basis comparison consists calculating the 
ratios for given conversion different values 
Various orders reaction are considered and the results present- 
Figures 27, 28, and for 90% conversion. 

The curves for the and models represent again 
limiting cases. Points lying the shaded areas correspond 
models which the zone perfect mixing located Figures 
20(a) and (b). 

inspection the figures shows that for orders 
superior unity, Model requires shorter residence time 
than Model MT, reach the same conversion. The reverse 
applies for reaction orders inferior unity. difference 

For the particular case values 6/6, have been 
tion, but MacMullin and Weber®® used the ratio 


For second-order reactions which curves show- 
ing the variation with the conversion have been presented 


Graphical Methods 


While design problems involving consecutive reactors can 
often solved analytically, long and tedious calculations may 
involved, the case Models and MT. order 
facilitate the solution such problems, graphical methods have 
been developed which involve, part, means computation 
whereby the overall conversion for number reactors 
series obtained from the individual conversions. 

Referring Figure 31(a), one can imagine 
reactor being made reacting zones series, each zone 
acting either tubular CST reactor. could also 
Figure 31(b). 

Since the one zone reactor becomes the feed 
the next, stepwise calculations are involved. 


for each the zones. Then, introducing Equation (4) one 
obtains: 

— Xe = 


Since combining the above expressions gives: 


But the overall conversion for the system, defined 
already, 


1—x = 


comparison with Equation (70), the overall conversion 
thus related the individual conversions the expression: 


1 —x=(l1 — x)(1 | 


may noted that this expression general because the 
particular conversion for any given zone reactor involves only 
the feed and effluent concentrations, irrespective the order 
reaction, the size the reacting zone even the degree mix- 
ing. 

interesting apply Equation (71) the particular case 
first-order reactions, for system consisting identical 
CST reactors series, the residence time each being 
The overall conversion for the system given by: 


72) 
which equivalent the expressions derived 
and Weber®® and Eldridge and the residence 
time was for the whole system, the residence time each 
reactor would and Equation (72) would 


the value this expression increases, the value 
becomes equal that Equation (74) for tubular reactor, 
indicating that, expected, infinite number CSTR 
series gives the same result tubular reactor the same total 
residence time. 

the application Equation (71), the conversion each 
reacting zone reactor determined graphically. This process 
facilitated means Figures 32, and 34, each which. 
presents together the curves for tubular and CST reactors. 
Similar curves, for one more orders reaction, have been 

For first-order reactions, Figure shows the variation 
obtained re-arranging Equations (5) and (12) respectively: 


Figure 33, for second-order reactions, plot Equations 
(6) and (13), while Figure 34, for third-order reactions, corres- 
ponds (7) and (14). 

general Equation (71) now applied Models 
becomes: 


which and represent the conversions the tubular and 
CST portions the reactor respectively. 


This equation used together with the following expression, 
which equivalent that relating partial residence times the 
overall for the reactor: 


where the residence time for the CSTR section and equal 
venience. 


Together, they make for easy solution the different 
required determine one the three quantities when 
two are known. Although details application vary slightly 
according the order reaction, especially for first-order 
reactions which the term Equation (77) out, 
the method will outlined way. Specific examples 
are worked out the Appendix. 
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Figure 31—(a) reactor with reacting zones series. (b) system 
comprising b different reactors in series. 


1.5 3.0 6.5 


Figure 33—Variation the conversion with CSTR and tubular 
reactors, for second-order reactions. 


Figure 35—Graphical construction for determining the conversion 
reactors. 


Determination the overall conversion when 
and are known. 


the analytical method necessary distinguish 
between models and MT. 


Model TM: 


The first step consists determining 
follows that: 


Since known, 


and 
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Figure 32—Variation the conversion with CSTR and tubular 
reactors, for first-order reactions. 
1.0 
0.6 
034 
0.0 
10 20 30 40 50 70 
Figure 34— Variation the conversion with CSTR and tubular 
reactors, for third-order reactions. 


Figure 36— Graphical construction for determining relative residence 


Depending the order reaction, read from the curve 
for tubular reactors, Figures 32, 34. noted 
that the abscissa these Figures contains the term 
refers the particular feed concentration the reacting zone 
reactor considered. 

this particular case, after conversion the first 
part the reactor, the feed the CST section will be: 


and thus the conversion will read abscissa 
the appropriate Figure. The values x,’ and 


are then introduced into Equation (76) give the overall 
conversicn 


1.0 
4 x 
0.6 
0.3 
0.0 
| 
— | 


Model MT: 


this case, read the CSTR curve the appropriate 

The concentration the feed the tubular section given 
by: 


and the overall conversion computed means Equation (76). 


Alternate Method: 


more simple method yet available for the special case 
reactors. consists graphical construction applicable 
either one Figures 32, 34. 

general way, one proceeds Figure 35, which can 
represent any one these three Figures. The conversion 
located the CSTR curve abscissa deter- 
mining point path then followed, parallel the 
portion the curve for tubular reactors, such that the abscissa 
the overall conversion for the reactor. 

similar graphical construction does not apply Model 
applies Model because can shown that 
segment equal segment CA. Indeed, the two sides the 
expression 


a 


represent segments and respectively. When verified 
mathematically, Equation (80) works out identity for 
other words, valid for all values between and 


Determination relative residence time, when 
and are known. 


Model TM: 


trial and error method involved which nevertheless 
yields results rapidly. Once value assumed, 6,’ intro- 
duced the term give corresponding value 
using the appropriate Figure 32, 34. Equation (76) 
then used compute the value which turn serves 
determine graphically the value of, When 

the sum 6,’ and thus obtained does not verify the 
value assumed originally for the procedure repeated until 
does. 

Since can read directly the curve for tubular reactors, 
using the appropriate Figure, the conversion considered, the 
ratio readily obtained. 


Model MT: 


procedure similar the above followed except that 
value first used, and then computed give 
Dividing this value xy)* yields 

the preceding case, the procedure repeated until the 
sum and equal the value assumed for 


Alternate Method: 


graphical construction again available for the particular 
case reactors. 

Figure represents curves for tubular and CST reactors, 
for any order reaction, Figures 32, 34. 

value assumed and the resulting abscissa 
determines point for the given conversion Following path 
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parallel determines the value the abscissa 
which, unequal the original one, requires that new value 
chosen and the procedure repeated until the value 
obtained agrees with the given one. 

can shown that the two sides expression similar 
Equation (80), corresponding segments and BD, are 
identically equal for all values between and 


Determination when and are known. 
Model TM: 


trial and error procedure indicated. Assuming value 
the corresponding conversion calculated according the 
method outlined already for determining the overall conversion 
when and are known. Should the result differ from the 
value given, the procedure repeated until agreement reached. 


Model MT: 


Although one can proceed with model TM, alternate 
graphical construction more suitable. 

Referring Figure 36, the values and determine point 
Following path parallel CD, point determines the 


Partial Mixing and Short-Circuiting 


The concept short-circuiting has been used differently 
statistical -passing channelling system. Short-circuiting 
reacting system refers, the present work, portion 
the feed which passes directly the outlet without having 
reacted been retained 

Referring Figure 12, the general equation for partial 
mixing and short-circuiting has been presented Cholette and 
Cloutier their Equation (8), which reads: 


Cr = Cu 
The expressions for the reaction rates and are: 


The concentration the zone perfect mixing, derived 
from material balance, given the following expression: 
Cy — (1 — 


(83) 
n 


Introducing Equations (11), (82) and (83) into Equation (81) 
the general expression reduces the following the case one 
reactant: 


or, terms conversion, 


XM 


may observed Equation (84) that the last term drops 
out when short-circuiting absent, whatever the order reac- 
tion. Furthermore, Equation (85), when and are 
equal unity. 

further observation shows that the term 
Equation (85) always positive always smaller than 
equal for the limiting case complete conversion. 
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Figure 37—Variation the conversion with reactors with partial 
mixing and short-circuit, for zero-order reactions, at residence times 
such that the conversion in a CSTR is xm, as indicated on curves. 
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Figure the conversion with and effect the short- 


circuit, reactors with partial mixing and short-circuit, 0.5, 
a residence time such that the conversion in a CSTR would be 90%. 


0.60 1.0 


— ne of x with m, in reactors with partial mixing and short-circuit, for first-order reactions, at residence times such that the 
CSTR would indicated fer each family curves. The effect the short-circuit, within family curves, characterized 
by the values of n. 


0.0 0.3 0.6 


Figure 40— Variation of the conversion with m and effect of the short- 
circuit, in reactors with partial mixing and short-circuit, fora = 1.5, at 
a residence time such that x4 = 90%. 


the latter case, and from Equation (83), 
given Equation (4), one obtains for the limiting case: 

Cr — (1 — n)Cr 
= 

While Fquations (84) and (85) are the general solutions for 
the model Figure 12, many particular solutions have been 
worked out, depending the order reaction and the magni- 
tude the short-circuit. For purposes comparison, the results 
are presented under two headings. 
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Variation the conversion with for given 
order reaction, different values the short- 
circuit. 


General Equation (85), applied different orders reaction, 
gives the following results: 
For such reactions, the short-circuit effect cancels 


The conversion thus dependent the volume 
mixing only. Equation (86), Figure 37, gives 
straight lines slopes all values 


(87) 


The curves Figure show that the short-circuit increases 
and the effective volume mixing decreases, the conversion 
departs more and more from the theoretical value, 


1.0 1.0 
0.9 
1.0 


0.20 


Figure 41—Variation of x with m, in reactors with partial mixing and short-circuit, for second-order reactions, at residence times such that 
the conversion xy in a CSTR would be as indicated for each family of curves. The effect of the short-circuit, within a family of curves, is 
characterized by the values of n. 
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Figure 42—Variation of the conversion with m and effect of the short- 
circuit, in reactors with partial mixing and short-circuit, fora = 3, ata 
residence time such that 


Families curves for first-order reactions are shown 
Figure 39. They are plotted different values residence 
time corresponding given theoretical conversions xy. 
Values ranging 0.4 are presented show the 
dominant effect the short-circuit higher conversions 


1.5 
Xu 
plot which shown Figure 40. 
a=2 
2 


Families curves derived from Equation (90) are presented 
Figure 41. for first-order reactions, the effect becomes 
dominant the theoretical conversion increases. lower 
values xy, the influence becomes more appreciable. 

When two reactants undergo second-order reaction and 
are present unequivalent molar amount, the following will 


apply: 


1.0 


0.6 


0.4 


Figure 43—Variation of the conversion with m, in reactors with partial 

mixing and short-circuit, for different orders reaction, residence 

times such that 0.9. Each family curves corresponds given 
value of the short-circuit. 


which reverts Equation (90) when 


XM 


plot which shown Figure for theoretical conversion 
Xu of 0.9. 


Variation the conversion with given 
short-circuit, for different orders reactions. 


Figure shows how the degree completicn reaction, 
conversion, affected the three variables and 

The presence short-circuit generally the determining 
factor affecting conversion, except low values 
certain value for given short-circuit, the mixing level has 
more and more pronounced effect the order reaction 
decreases; high values are reached, the conversion 
becomes approximately the same for all orders reaction. 

Generally, observed that for given values and the 
conversion increases with increasing orders reaction all 
values smaller than The reverse true for values 
larger than When the conversion the same for all 


(91) orders reaction, indicated the curves one family 
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CONCLUSIONS 


The present study should help throw light the perfor- 
mance chemical reactors general because, pointed out 
Smith correlations were available yet the chemical 
engineering literature for partially mixed reactors. Although 
rather simple models have been investigated, there should 
follow better understanding the relative importance differ- 
ent factors affecting the behavior actual reactors. 

further investigations are carried out, other 
will considered and addition more factors may have 
taken into account such actual mixing patterns, turbulence, 
diffusion and mass transfer from one zone another. 

While this study restricted homogeneous, irreversible 
and isothermal reactions, other interesting results are anticipated 
when reversible, consecutive, simultaneous reactions are 
carried out isothermally adiabatically, models such those 
defined the present work. 

evident that the scope chemical kinetics yet 
unlimited and that great number proinising research topics 
should investigated order secure sound information 
the phenomena intimately associated with chemical 

Although reactors must designed and built even when 
data rigorous design procedures are not available the ever 
increasing contributions the subject add the scientific and 
mathematical background reactor design. 
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Nomenclature 


refers reactant 

refers reactant 

volumetric concentration, 

effluent concentration. 

feed concentration the reactor. 


feed concentration the second zone reactor, 
models MT. 


reactants leave one reacting zone and enter another. 
theoretical concentration the effluent mixing was 


perfect. 

theoretical concentration the effluent the reactor 
was tubular. 


velocity rate constant 

fraction the total volume which perfectly mixed. 
fraction the feed entering the zone perfect mixing. 
volumetric rate flow, 

reaction rate, 


reaction rate, CSTR. 

reactor. 

reaction rate for the well-mixed zone reactor. 


fraction the total volume undergoing piston flow (special 
case: Figure 20). 
volume, 


theoretical conversion ina CSTR. 

theoretical conversion tubular reactor. 

reactor. 

reactor. 


order reaction. 


residence time, hr. 

model. 

model. 


identical series gase: Equation 
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APPENDIX 


few examples are worked out illustrate much 
possible the various applications the graphical methods des- 
cribed, changing not only the different parameters but also the 
form solution. 

(1) Find the conversion for Model MT, given: 


(0.3) (2.3) 0.69 
and from Figure 32, 0.41 
(0.7) (2.3) 1.61 
From Figure 32, 0.8 


From Equation (76): 
0.88 


Since for 2.3, Figure gives: 0.9, the actual 
conversion 0.88 found above can checked Figures 

(2) Find using graphical construction, for Models 

Point the general construction Figure 35, found 
have the following coordinates, when using the proper curves, 
Figures 32: (2.3, 0.88). 

Following path parallel DC, the abscissa point 

The value then 0.69/2.3 0.3 

(3) Find for Model TM, given: 


9.0, 


from Figure 33, 0.47 

From Figure 33, 0.62 
Equation (76) gives: 


the reactor was tubular, value 9.0 would give 
0.9; the above conversion, 0.8, can checked 
Figures 14, 24. 

(4) Find for Model TM, given: 0.5, 


From Figure 33, 0.2 
(0.83 
0.83) 


Figure gives then: 0.375 


1.22 0.2) 0.98 


Equation (76): 


which checks the given conversion. Had been different, 
new value would have assumed and the proce- 
dure repeated until the conversion obtained would correspond 
the given one. 

Now, the total value is: 


0.25 1.22 1.47 


The required ratio then 1.47/1 1.47 which checks 
the value given Figure 16. 

(5) Find 6/6,, using graphical construction, for Model 
MT, given: 


Assume 

Using Figure and applying the construction Figure 
determines the value 58.4/99 0.59, which corresponds 
the value given originally. had been different, other values 
would have been assumed, until the resulting value 
agreed with that given. 

obtain conversion 0.9 tubular reactor, Figure 
shows that 49.5 

The ratio then: 99/49.5 which value can 
checked Figures and 30, 

(6) Find for reactor, given: 


Thus, (2) (49.5) 99. 


Figure gives 0.76 and Equation (76) yields: 
0.9 


0.76 


= 1 


40.6 
= 


The total value is: 58.4 40.6 and checks 
the value derived first from the given ratio 

not the same, another value would have 
assumed and the procedure repeated until agreement was reached. 

The value obtained from the ratio: 


58.4 


Figures and verify this result. 
graphical construction, example would give the 
same result without involving trial and error calculations. 
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The Incorporation Fission Products 


into Glass for Disposal’ 


BANCROFT? 


There need economical and safe method 
disposing radioactive waste produced nuclear 
reactors. method for the incorporation fission 
products into glass being studied Chalk River. 
Glass containing curies fission products 
per kilogram has been made during the investigation 
the process. This radioactive glass has been used 
stability studies the laboratory and the field. 


Nepheline syenite rock containing Na, and 
silicates) fluxed with 30% weight lime, 
used the base for the glass. nitric acid solution 
the fission products added pelletized mix- 
ture the solids. When the mixture heated nitric 
acid evaporated near 110°C. The nitrates are de- 
composed the temperature range and the 
mixture melted form glass maximum 
temperature 1350°. Two the fission products, 
ruthenium and cesium, are volatilized varying 
degrees from the mixture during heating and create 
major contamination and containment problem. 


The method processing and results obtained 
are discussed. The volatilized elements have been 
removed from the off-gas stream method which 
permits their incorporation into glass. Scale-up 
the process treat fission product wastes from 
power reactors should present great difficulty 
other than design equipment for the remote pro- 
cessing large quantities radioactive material 
high temperature. 


fission products generated nuclear reactors are 
biological hazard because they are radioactive. This 
radioactivity eliminated only natural decay, which 
occurs times ranging from seconds many years for 
the different nuclides. Two the longer-lived nuclides, 
strontium-90 and cesium-137 represent hazard for 
hundreds years. 

Spent fuel from nuclear reactor may processed 
chemical extraction methods separate the uranium 
and plutonium from the fission products. The fission 
products are concentrated nitric acid solution and 
this solution stored waste metal tanks. this 
method storage were continued for the whole the 
hazardous lifetime the fission products would 
necessary transfer the solution periodically new 
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tanks because for long periods, the tanks cannot guar- 
anteed against leakage. The cost this storage high. 


the future, when many nuclear reactors are used 
generate electric power, the rate production 
fission products will very much higher than now. 
The cost treating the radioactive wastes for disposal 
will not represent large percentage the over-all cost 
the electric power, but may appreciable percentage 
the cost reprocessing fuel and must kept low 
possible; hence the need for economical and safe 
disposal method for fission products. safe method may 
defined one that does not allow significant 
fission products escape man’s environment for 
many hundreds years. 


There are many methods disposal being considered 
groups different countries. Some these are listed 
follows: 


solution 
(1) metal tanks 
(2) deep underground strata 
(3) underground salt formations 


solids 
(4) solution converted calcined solid 
(5) fission products incorporated into fired clay 
(6) fission products incorporated into glass. 


Work being done Chalk River develop 
method incorporating fission products into glass. Some 
the desirable features glass medium for disposal 
are follows: 


greater variety burial locations are available 
than for methods (2) and (3). 

There less risk widespread dispersal radio- 
activity from solid disposal than there from 
liquid disposal. 

The rate release fission products water 
attack lower from glass than from the 
calcined alumina (method 4); there not enough 
information available make comparison with 
fired clay. The rate dissolution water 
important factor when considering solids, because 
buried the ground, they may contact with 
ground water. they products 
water unsafe rate they would have 
stored water-proof containers, such metal cans 
concrete vaults, and this extra protection would 
add the cost the disposal. expected that 
glass can made such quality that can 
buried without added protection. 


Several years ago was decided that experimental 
equipment should handle one gallon batches 
fission product solution. The purpose building this 
equipment was twofold: 


prepare pieces glass containing sufficiently 
large quantities radioactivity allow reliable 
tests done the laboratory and the field 
the dissolution glass water. 


obtain information operating problems, such 
dusting, volatilization, handling and control. 


During 1956 and 1957 equipment development, construc- 
tion and testing were carried out, and May 1958 the 
first block highly active was produced. Since that 
time more than 2000 curies fission products have been 
incorporated into glass, most which being used 
leaching tests. 


Cost estimates were made for incorporating the fission 
products concentrated power-reactor wastes into glass 
the process described this paper. Processing costs 
were the order $0.05 per gram fission products 
per kilogram glass. These costs represented less 
than the cost the electric power produced, but 
the additional expenses disposing the glass were 
considered increase the over-all disposal cost 1-2% 
the power cost. 


DESCRIPTION EQUIPMENT 


Figure schematic flowsheet the process 
shown. The solid and liquid ingredients are mixed the 
crucibles which the glass melted. The crucibles 
are transferred from the mixing station the drying 
and denitrating furnace, which they are heated 
maximum temperature 900°C. After this heating the 
mixes are allowed cool over night, and then are placed 
the melting furnace for further heating, 1350° for 
one hour. The glass which formed cooled the 
furnace before removal for testing. The gases generated 
both these furnaces are drawn through gas-cleaning 
system where the volatilized radioactivity collected 
adsorbers, the nitric acid recovered packed column, 
and the remaining gases are passed through caustic 
scrubber and filter before being exhausted. 


Before description given the equipment used 
for this work should mentioned that, because these 
fission products are highly radioactive, they must 
handled with extreme care. Adequate shielding (four 
inches lead two feet concrete for this equipment) 
must provided protect personnel from excessive 
exposure radiation. All equipment used the cells 
must reliable and simple maintain, because repairs 


FISSION 
PRODUCT 
SOLUTION 


NEPHELINE 
SYENITE 


Ru Cs 
ADSORBERS 


DRYING 
FURNACE 
900°C 


can often made only remote methods. All handling 
systems must suited the use hand-operated 
manipulators and pneumatic, hydraulic and electrically 
operated mechanisms. Thus developing process for 
handling radioactive materials attention must given 
each operation. The simplest means are usually the best. 


Figure shows some the equipment the operating 
cell. This photograph was taken some time before the 
system was completed and not accurate all respects. 
guide size, the crucibles are in. diameter 
the top and inches high. the left the photograph 
the mixing station, box that contains the crucible 
while the solution and solids are added. The shroud 
lowered cover the crucible while the solution being 
added, and reduces the spread activity splashing. 
Liquid added through the pipe entering the front 
the box and pours into the middle the crucible. The 
pelletized solids are also fed the middle the crucible, 
through in. pipe which not shown. The box fitted 
with door moved hydraulically cover the opening, 
has water spray for washing, and ventilated 
separately from the rest the cell. 


Each furnace can accommodate two crucibles. The 
furnaces are home-made and are sections for easier 
disposal; they are heated eight silicon carbide heating 
elements. Electrical connections the elements are made 
wires through the ceiling and spring- 
loaded contacts under the furnaces. These connections 
are necessary allow the remote replacement ele- 
ments. The doors are moved electric hoists. Inside the 
kw. drying furnace, the left, stainless steel can, 
which serves contain the nitric acid fumes and vola- 
tilized activity. cover for the can fitted the inside 
the furnace door and can pressed against the can 
hydraulic piston. The melting furnace (30 kw.), 
the right, contains metal but are two shallow 
silicon carbide receptables which are large enough 
contain the glass from crucible cracks. 


the centre the photograph elevator which 
travels through hole the ceiling and used bring 
crucibles into the cell. Figures and show the handling 
tongs used move the crucibles about the cell. Once 
the glass has been prepared moved means 
dolly track into adjacent cell. This cell has pro- 
vision for storing glass, removing from the crucibles, 
leaching flowing static water, and removing 
from the cell for disposal. Figure shows crucible being 
placed into leaching tank. The tongs are moved through 
ball-and-socket joint the ceiling and have air-operated 


jaws. 
causric 
SCRUBBER 


FISSION 
PRODUCTS 
IN GLASS 


Figure 1—Flowsheet for the incorporation fission products into glass. 
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Figure 4—Handling tongs outside 4-inch thick lead shield- 
ing wall, showing window and ball-and-socket joint for 
tongs. 


DESCRIPTION PROCESS 
Mixing 
The composition nepheline syenite (N.S.), natur- 
ally occurring rock used the glass, given Table 


TABLE 
APPROXIMATE COMPOSITION NEPHELINE SYENITE 


Si0, 60% 10% 
24% 


others 


Nepheline syenite can melted the temperature range 
1300-1400°C., but extremely viscous and cannot 
made free bubbles convenient heating time. How- 
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Figure 3—Handling crucible inside cell with tongs. 


Figure 5—Placing crucible into leaching tank. 


ever, lime added the N.S. bubble-free melt can 
made heating for times short one hour 1350°. 
used, but the normal mix 85% N.S.-15% CaO. 

The waste solution which was used for the work 


described this paper had the composition indicated 
Table 


Figure 2—Operating cell equipment. 


TABLE 


ANALYSIS CONCENTRATED FISSION-PRODUCT SOLUTION 
February, 1959 


Acidity FINO; 

Specific gravity 1.39 

Total 8.6 10° counts per minute per millilitre 


227 mg/ml (Na, Fe, Ca, Al, Cr, Ni, Mg) 
100 mg/ml 


Total solids 
Non-volatile solids 


When the solution and solids are mixed together the 
nitric acid decomposes some the N.S. and forms 
silica gel. The gelling time for mixture this solution 
and pure N.S. about seconds. 15% lime mixed 
with the N.S. before the solution added the gelling 
time somewhat longer, depending upon the method 
mixing with the acid solution. Some heat generated 
the reaction and this can cause the temperature the 
mixture increase the boiling point HNO,. 


During the early development the process using 
solutions which were not radioactive, rotating stirrer 
was used mix the solids and liquid. For some mixtures 
uniform blends were obtained stirring for one minute 
and the stirrer could removed before the mixture 

elled. However, with some compositions the mixture 
gelled very short time and the stirrer was either 
trapped had large fraction the mixture adhering 
removal. stirrer was therefore unsuitable for 
remote operation and the powdered nepheline syenite and 
lime are now mixed together and pelletized rolling 
rotating drum. The pellets are screened and the +20 
mesh fraction used. this form the solids can fed 
from vibrating feeder, and mixed with the solution 
without stirring agitation. The pellets have enough 
free space between them allow the solution soak 
through layer several inches thick before gelling blocks 


the channels. 


For convenience operation the mixes are prepared 
feeding layer pellets several inches deep into 
crucible and then adding enough solution react. Alter- 
nate additions solids and liquid result layered 
mixture which has dry pellets the top layer. 


The silica gel formed mixture generally quite 
firm, and after being dehydrated heating very 
porous. This texture the gel reduces the possibility 
some the mixture being “bumped” from the crucible 
during heating. addition, the porous structure the 
dried mixture provides channels through which gases leave 
the mix. These channels are tortuous and are probably 
efficient de-entraining particles from the gases flowing 
them. The entrainment losses from mixture are less 
than 0.1%. 


Drying and Denitrating 

When mixture solution and solids heated many 
reactions take place. temperatures 100-110°C. the 
water and nitric acid boil off. The silica gel dehydrated 
the range 100-150°. 135° aluminum nitrate de- 
composes release oxides nitrogen. higher temper- 
atures the nitrates iron, sodium, potassium and calcium 
decompose, the latter about 600°. Since the temperature 
the mixture not uniform during heating, probable 
that many these reactions occur The 
maximum temperature the first phase the heating 
900°, which high enough ensure complete denitra- 


tion, but does not cause rapid deterioration the stainless 
steel can inside the furnace. 

From the furnace the vapors and gases released during 
heating and the air entering through cracks are drawn 
through gas-cleaning system. The 
activity collected granular adsorbent and the nitric 
acid and water are condensed column packed with 
Raschig rings. The oxides nitrogen are scrubbed the 
same column with nitric acid solution convert them 
nitric acid. After passing through the condenser- 
recombiner the gases are bubbled through solution 
sodium hydroxide remove entrained acid. They are 
then passed through filter and exhausted the building 
ventilation system. 


Melting 

The second phase the heating done furnace 
which can heated 1400°C. The gases from this 
furnace are also passed through the packed column, which, 
this stage, operated gas scrubber rather than 
condenser—recombiner because most the volume air 
drawn through cracks the furnace. 


After the mixes have been heated 1350° and kept 
that temperature for one hour, they are cooled slowly 
inside the furnace anneal the glass. The temperature 
after hours about 300° and after another about 
room temperature. After this time the glass removed 
for inspection and testing. 

photograph 2-kg. piece 85% N.S.—15% CaO 
glass shown Figure This glass clear and almost 
colorless, but that made with fission product solution 
dark brown, colored mainly iron and chromium cor- 
rosion products. The upper surface free from crust and 
foam, but there are few surface bubbles evident. Parts 
the curved surface were chipped away when the 
crucible was broken from the glass and some the 
crucible clay was left adhering the glass. The physical 
features this specimen, except for color, are typical 
the glass produced during this work. 


Volatilization Losses 

Two the major fission products are volatile under 
conditions that exist during the heating mix. Ruthen- 
ium and cesium volatilize and escape varying degrees 
from the mixes. This presents the most troublesome con- 
tamination problem the process. 

mixture containing free nitric acid most the 
ruthenium which volatilizes escapes temperatures 


Figure 6—A 5-inch diameter hemisphere 85% N.S.— 
15% CaO glass. 
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below 900°. the acid neutralized with lime then 
the occur mostly above 900°. The losses from 
both the acid and neutral mixes are about the same for 
complete heating cycle. the present equipment the 
volatilized activity can best contained the lower 
temperature range and the mixes are made contain 
free nitric acid. 


indicate that most the cesium that volatilizes leaves the 
mix the temperature range 400-700°, but losses continue 
occur the maximum heating temperature 1350°. 


Ruthenium and Cesium Collection 


course, would desirable keep the ruthenium 
and cesium the mix, but suitable method doing 
this has been found. Typical volatilization losses from 
kg. mix heated 1350° are 50% the ruthenium and 

the cesium. the solution used, ruthenium repre- 
sented about the total radioactivity, and the cesium 
was about 40%, the loss from mix was about 
the total activity. Two per cent the total represents 
appreciable amount radioactive material and this 
could not allowed move uncontrolled through the 
equipment. was therefore necessary collect this 
volatilized activity one location and such form that 
could handled conveniently for disposal. 


was noticed laboratory experiments that ruthen- 
ium was concentrated from the vapor the stainless steel 
surfaces the apparatus. was believed that iron, 
oxides iron, had great capacity for retaining ruthen- 
its oxides. Experiments were done test iron 
oxide adsorber ruthenium. provide high 
surface area, porous firebrick that had been crushed 
size +20 mesh was soaked solution ferric 
nitrate. The firebrick was dried and the nitrate decom- 
posed heating, leaving coating ferric oxide 
the granules. Ruthenium volatilized from mix, along 
with nitric acid vapor and air, was passed through bed 
the granular adsorbing material 4-in. deep. About 98% 
the volatilized ruthenium remained the adsorber. 
The same collecting efficiency was observed for adsorber 
temperatures 1000° and for superficial gas velocities 
one foot per second. This type adsorber has 
been used successfully the vapor outlets the two 
furnaces. 


Cesium was not adsorbed iron well was 
ruthenium, but did collect any surface that was cool 
enough allow condense. bed granulated fire- 
brick coated with iron oxide could, therefore, used 
collect cesium part were below 400°C. 


future plant the adsorbing material used collect 
most the activity volatilized from mixes will contain 
only small percentage the through-put 
product activity. However, this amount activity will 
appreciable and should not considered disposal 
separate from the glass. The present intention incor- 
porate the adsorbing material into the glass. This addition 
will change the glass composition and may reduce the 
resistance attack water. This problem under 
investigation. 


LEACHING TESTS GLASS 


One the advantages glass medium for dis- 
posing fission products its extremely low rate 
dissolution water. Therefore, important that this 
property the glass measured and related the 
factors which effect it, such composition, conditions 
preparation, and conditions exposure water. 
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For several years, leaching tests have been conducted 
small specimens glass (less than grams) con- 
taining radioactive nuclides. Most these tests were 
done immersing the glass distilled water polythene 
envelopes. Each day the specimens were transferred 
fresh water new envelopes. The amount activity 
released from the glass was measured counting 
sample the water, and the leaching rate was expressed 
grams Although was known that 
all the components the glass did not leach the same 
rate, this method expressing apparent leaching rate 
was simple and adequate for comparing glasses different 
compositions. The lowest leaching rate measured these 
tests was for nepheline syenite- 
15% lime glass containing Ru-106 the tracer nuclide. 
For glasses other compositions the rates were high 

The major difficulties doing leaching tests were 
found sampling and analysis. Since the rate 
attack glass was very low, the concentration leached 
material the water was low 0.005 ppm. result 
this low concentration, inconsistencies were found 
testing the small samples and was thought necessary 
improve the method. 


The tests now being done flowing water systems 

are improved the following ways: 

The glass has higher specific activity. 

The ratio surface area glass water volume 
higher. 

The glass specimens are much larger and have 
surface areas 500 times those the smaller 
specimens. This reduces the possibility minor 
disturbances affecting the tests appreciable 
extent. obtain large surface areas other workers 
have used powdered glass for leaching tests. Two 
the major difficulties for this technique which 
make undesirable are (a) determining the surface 
area, and (b) separating the entrained solids from 
the leaching solutions. 


Figure plot leaching rate against time for 
two separate tests groups four 2-kg. blocks. This 
shows initially rather high rate which decreased 
These rates are typical nepheline syenite glass con- 
taining 15% lime and 1.5% solids from the fission product 
solution. The amount fission products released from 
glass leaching these rates small that direct burial 
the ground might acceptable. 


RATE 


DAYS FROM START LEACHING 


Figure Leaching rate glass. Nepheline syenite 
15% lime. 
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Currently glasses different compositions are being 
leached. Tests are also being made the laboratory 
measure the retention leached nuclides 
Twenty-five 2-kg. pieces glass containing about 300 
curies fission products have been buried the field. 
Samples the ground water nearby are taken period- 
ically, but not enough data are yet available estimate 
the rate release fission products from this glass. 


CONCLUSION 


Glass, stable solid, desirable medium for the 
disposal products. the glass exposed 
water while buried the ground should release 
fission products only extremely low rate. The pro- 
cess being studied Chalk River for the incorporation 
fission products into glass appears suited large 
scale operation. Development studies will continued. 
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Simplified Method for Determination 
Solid Diffusion Coefficient with 


Non-Linear Adsorption 


CHI TIEN? 


simplified method proposed which enables 
the determination solid diffusion coefficient for 
systems with non-linear adsorption isotherm. The 
method involves the approximation the concen- 
tration history the liquid phase with polynominal 
functions and consequently the match the experi- 
mental results and calculated values. Experimental 
work has given rather satisfactory results. 


the study rate processes heterogeneous systems 
involving solid and fluid, like adsorption ion 
exchange operations, the knowledge the diffusion co- 
efficient the solute into the solid phase primary 
importance. This especially true the solid diffusion 
the dominant factor the overall mass transfer mechanism. 

One the several common methods which are used 
for measuring the diffusion coefficient solid the 
sorption from constant and limited volume well- 
stirred fluid. this method the progress sorption can 
easily followed. the other hand, solution the 
diffusion equation for the appropriate boundary condi- 
tions, leads theoretical expression for the progress 
sorption terms parameter D.t. Comparisons the 
theoretical curve and actual curve would yield the numeri- 
cal value the diffusion coefficient. The detailed 
description this method has been reported several 
publications 

One restriction imposed this method the distri- 
bution coefficient the solute the fluid and solid phases 
has constant or, other words, the adsorption 
equilibrium relationship should linear. This necessi- 
tated the fact that only the case constant distri- 
bution coefficient can the solution the corresponding 
diffusion equation obtained. Unfortunately larger 
number systems exhibit non-linear behavior. has 
been suggested that case non-linear systems 
differential technique could used. However, this 
method, just like any other differential method, requires 
precision set-up and often gives unreliable results. 


The object this investigation present method 
which enables the determination solid diffusion coeffi- 
cient for non-linear system the utilization the 
information sorption from finite and constant volume 
well-stirred fluid. This accomplished substituting 
dependent variable independent boundary condi- 
tion. The detailed description given the next section. 
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Mathematical Development 

The diffusion solute from finite and constant 
volume well-stirred fluid into initially empty spherical 
particles can described follows: 


(4) 


The assumptions involved are: (a) the spherical par- 
ticle assumed quasi-homogeneous structure, 
(b) the diffusion constant, assumed constant, 
and (c) the desorption process insignificant. The word, 
well-stirred, means that the concentration fluid 
phase uniform any time and, furthermore, because 
the rigorous agitation, the resistance the diffusion 
process due the liquid film negligible. This also 
implies that the solute the surface the particle 
with that the fluid phase. 


Solutions Equation (1) with boundary condi- 
tions given Equations (2), (3) and (4) have been 
given Carslaw and Wagner and 
Wilson for the case constant distribution coefficient. 
For the case where function concentration, the 
resulting integral equation becomes non-linear and the 
solution becomes not obtainable. 

can seen from Equations (1), (2), (3) and (4), 
complete would give the expression the 
solid phase concentration, function time and 
position (inside the spherical particle) and expression 
the fluid phase concentration, function time. 
The relationship between and given Equation 
(4) and, furthermore, neither them independent. 

This problem, however, can greatly simplified 
the fluid concentration known. that case, with the 
assumption equilibrium between the fluid phase and 
the surface solid particle, can obtained con- 
junction with the particular absorption isotherm. Further- 
more can approximated function with poly- 
nominal expression degree 

m 

should noted the coefficients, can only 
determined from direct experimental observation. 

With given Equation (5), solution the ori- 
ginal diffusion equation found be: 
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Let the amount solute adsorped per unit mass solid, given the following expression: 
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Substitute Equations (9) into (10) and after rearrangement have 
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With the surface concentration known, can 
plotted against time with a?/D.t parameter. This would 
result series curves. the other hand, from the 
experimental results, can plotted against time. 
This experimentally obtained curve will intersect the 
computed curves number points. Each these 
points specified particular values and a?/D.t. 
This turn gives the numerical value Furthermore, 
consistant value would prove the validity the 
assumption that the diffusion coefficient constant. 
illustrative example will given the next section 
demonstrate the detailed procedure. 


Illustrative Example 


The adsorption oxalic acid from glycol solution 
with Permutit SKB has been studied For particular 
run, the following results were obtained: 


500 cc. glycol solution with 0.0512 millequivalents 
oxalic acid per cu. cm., 2.7 grams (dry basis) Permutit 
SKB were added and vigorously agitated eliminate the 
resistance the liquid film. The resin when completely wet 
was 0.049 inches diameter. The concentration the liquid 
various times was from titration. 


1000 0.665 4.64 0.700 
4500 0.524 


Values were obtained with the information 
adsorption isotherm the system. q/q, were found 
material balance. 

The surface concentration q,, function time, 
approximated polynominal second degree. 
The equation found be: 


The comparison between the actual value and 
that given this approximation presented Table 
Generally speaking, values given Equation 
(15) are lower than the actual value for the value 
time and are higher for higher values time. 
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TABLE 


Comparison between actual surface concentrations and values 
given approximation. 


Surface Concentration, 


Actual Approximate 
300 4.62 
600 4.55 4.5 
1000 4.47 4.455 
1800 4.33 4.335 
3000 4.26 4.31 
4500 4.24 4.27 


With given (15), (11) 
reduces to: 


dav 1 6 


-D 
> P; ti 
i=0 i=) 
-D 
P; ti 
i=J) 
4.24 


Numerical values q/q, corresponding different time 
and preassigned value D.t/a? were computed from 
Equation (17) and graphically presented Figure 
the same chart, the curve relating q/q, and time and 
observed from the experimental work was superimposed. 
This latter curve intersects the calculated values q/q, 
the following points: 


500 1000 1500 
Time seconds 


2000 


Figure 1—Analysis experimental data for determination 
diffusional coefficient. 


i=0 | 


Time Sec. a?/D.t 
482 
0.0245 in. 


this particular case, was found increasing 
with time, however, other cases, the opposite trend 
was observed. Since definite trend exists, simple 
arithmetic average was taken. felt that least part 
the deviation due the inaccurate approximation 


Conclusions 

simplified method proposed which enables the de- 
termination the solid diffusion coefficient for systems 
with non-linear adsorption isotherm. This method requires 
very simple experimental set-up and the results from 
limited work seem rather satisfactory. 


Nomenclature 


Radius spherical particle 

Concentration fluid phase 

Initial concentration the fluid phase 

Diffusion 

Product the surface area solid particle and its 
density 

Coefficients the i-th term the polynomial expression 
for surface concentration 

solid phase, based per unit mass 
soli 

Amount solute adsorped, per unit mass solid 

Surface-concentration, based the per unit mass 

Radial distance from the centre the spherical particle 

Time 

Volume total fluid 

Distribution coefficient 
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Nuclear Grade Uranium Tetrafluoride 


the Moving Bed 


HUESTON? 


After two years research and development 
the pilot plant stages, plant for the production 
nuclear grade uranium tetrafluoride the Moving 
Bed Process has been designed, built, and successfully 
operated Eldorado’s Port Hope Refinery. This 
the only production facility the world using the 
moving bed technique for the production 


exceptionally high grade uranium tetrafluo- 
ride, low unconverted uranium oxides and uranyl 
fluoride can obtained. large percentage 
the production exceeds 99% uranium tetrafluoride. 
Total metal contamination less than 150 ppm. 


SUPERHEATED 


cess from uranium UO, produced 
Port Hope the Solvent Extraction Process (See 
Figure 1). 

The first step the hydration the UO, and the 
casting the resultant into pellets; next comes 
the reduction the hydrate UO, moving bed 
reactor countercurrent flow hydrogen; and finally 
there the hydrofluorination UO, with counter- 
current flow gaseous hydrofluoric acid. 


1Manuscript received June 9; accepted October 22, 1959. 
2Development Engineer, Eldorado Mining and Refining Limited, Port 
Hope, Ont. 
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Figure 1—Production uranium tetrafluoride moving bed process flow sheet. 
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Previous Work 


Original laboratory studies the feasibility this 
process are described Moore 

Laboratory bench scale studies were carried out 
Brunswick, New Jersey. 

great deal work the evaluation the moving 
bed process was performed the U.S.A.E.C. plant run 
National Lead Company, Fernald, Ohio.. summary 
this work was given Arnold the recent 
Geneva ‘Atoms for Peace’ Conference. The project 
Fernald was shelved 1957. 


pilot plant with two 8-in. diameter reactors was 
constructed the Port Hope Refinery the fall 
1956 and operations started December the same year. 
Melvanin discussed the layout and operating conditions 
the plant. 


took about months overcome the production 
problems; and the Spring 1958 was decided 
construct production scale plant. very tight schedule 
was laid out, and Catalytic Construction Company, the 
construccor the pilot plant, was retained design and 
build the new plant conjunction with Eldorado per- 
sonnel. The new plant started December 11, 1958 and 
produced specification product before the end the year. 


Theoretical Discussion 


not the intention here discuss the general theory 
behind moving bed reactors, this has been adequately 
covered elsewhere 


The first step the process the formation 
hydrate pellets according the following reaction* 


The reason for making the pellets two-fold. Firstly, 
feed discrete form must fed the reactors 
avoid packing. Secondly, research investigations have 
shown that the hydrate more reactive feed material 
than the UO,. One theory that the hydrate pellet has 
its moisture removed the first stage treatment and 
left quite porous, therefore offering large surface 
area for gaseous contact subsequent reaction: 


The third step may accomplished one two 
ways. The first direct hydrogen reduction: 


The second two stage process—a thermal reduction 
followed the hydrogen reduction: 


+ 3U0, + 2H:0 + Heat (5) 


There very little control over which way the reduction 
takes place. Gas consumption data indicates finally 
accomplished partly shown both reactions. Some 
typical figures are given Table 
The fourth step hydrofluorination. The main and 
desirable reaction the hydrofluorination reactor is: 
undesirable side reaction, caused unconverted UO, 
U,O, forming urany] fluoride also takes place: 


hydrates uranium trioxide are discussed Bright and Jongejan 
(8). questionable whether the. only hydrate present. 


TABLE 
Gas CONSUMPTION DATA FOR THE CONVERSION UO; 


Percent theoretical 
basis UO; UO, 


Flow Dissociated 


Product Rate Ammonia Cu. Ft. 


Ibs. Per. Hr. per. H.* Reduction 
175 285 
220 360 93.0 
270 470 98.2 


*Reduced N.T.P. 


Because these side reactions, extreme care must 
taken complete the reduction phase, and protect 
adequately the hydrofluorination reaction from the air. 
This done blanketing the reactor with inert gas. 


Eldorado Development 


One the greatest problems the pilot phase was 
the development equipment which would make pellets 
sufficient reactive qualities produce high grade UF, 
from the gas-solids reactions follow. 


Previous work the U.S. Henline indicated that 
cast pellets appeared highest reactivity with 
regard the conversion Preliminary work 
Port Hope indicated that this was so, much the 
commercial pellet machines that were tried produced 
unsatisfactory pellets. every case, the unequal hardness 
pellets due concentric layering and/or too high 
density caused poor conversion from UO, UF, 
inhibiting the diffusion the anhydrous hydrofluoric “acid 
throughout the pellets. Eventually, became necessary 
abandon the commercially available machines and develop 
separate techniques and equipment. 


The preliminary design the reduction reactor 
proved quite good and only minor changes were 
required make this process trouble-free. 


The major problem the hydrofluorination reactor 
was the maintenance continuous flow through the 
equipment. During the reaction from UO, UF, 30% 
expansion takes place and tends cause bridging the 
reaction rate becomes too rapid any one time, any 
portion the reactor. There also tendency for the 
pellets soften during the conversion period. The placing 
vibrators key locations the reactor walls and the 
development proper operating techniques greatly re- 
duced the plugging tendencies the reactor. 


DESCRIPTION EQUIPMENT AND PROCESS 


Pelleting Equipment 


Uranium trioxide powder continuously fed 
adjustable screw feeder into mixing vessel into which 
controlled volume water also fed. The resultant 
flows continuously onto rubber casting belt, 
which approximately ft. long and runs ft. per 
minute through controlled climate cabinet maintained 
temperature 180°F. and relative humidity 
+90%. the end this belt the pellets are knocked out 
their cavities rotating mechanical device. They 


discharge onto curing belt moving approximately 2/3 


ft. per minute opposite the casting belt. 
this operation, temperature 70-80°F. maintained, 
and relative humidity set facilitate removal excess 
moisture from the pellets. (See Figure 2). 
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Figure Pellet casting and curing equipment. Note: 
Screw feeder for oxide additions and flowrator for the 
water. Part dust collection duct work may seen 


the far end. 


The whole pelleting operation enclosed dust- 
tight cabinet with dust facilities prevent 
contamination the plant air with uranium dust. 


The pellets are pyramidal diamond shape, approxi- 
mately 5/8 in. long, 3/8 in. wide and 3/16 in. deep the 
base. The unit has capacity 1000 hydrated 
pellets per hour. This sufficient maintain 24-hour, 
seven-day design capacity operation the reactors 
operating 16-hour, five-day week. 


Reduction Reactor 

The reduction reactor 14-in. internal diameter 
vessel, approximately ft. high. has gas disengaging 
section located ft. from the top the reactor for the 
release off-gases from the resulting reactions. (See 
Figure 3). Design capacity one mole per hour UO, 
(270 Superheated steam 600°F. and in. W.C. 
pressure admitted concurrently the top the reactor 
along with the pellets. This dehydrates the pellets for the 
gas-solids reaction follow. Because the presence 
small quantities nitric acid, the top ft. the reactor 
type 347 stainless steel. The remainder the reactor 
309. 


The actual reduction from UO, UO, carried out 
the lower portion the reactor 1000-1100°F. The 
reducing media dissociated “cracked” ammonia (3H, 
which generated passing anhydrous ammonia 
over nickel catalyst 1700°F. Temperature the 
reaction controlled the introduction blanketing 
gas which analyzes +90% nitrogen, carbon dioxide 
and small quantities hydrogen and carbon monoxide. 
This gas produced the controlled combustion 
natural gas with air furnace, with the product gases 
being collected and compressed 150 for use. 

The reactor heated three jackets, each 
with maximum output kw. The power the 
jackets set maintain stipulated bed temperatures 
throughout the reactor. unique switching arrangement 
allows the choice six levels heat maintain condi- 
tions and prevent heavy on-off demand for power and 
consequent cyclic thermal conditions the reactor wall. 
Table gives some typical data with regard operation 
the reactor. 
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Figure General view reactors. Hydrofluorination 
reactor the right and reduction the left. 


OPERATING FLOWS AND PRESSURES THE REDUCTION REACTOR 
VARIOUS THROUGHPUTS 


Gas velocity in. W.C. 
Product rate| rate ft. per sec.* reactor 
lb. per hr. ft. per hr. (void space) disengaging section 
125 0.70 0.50 
0.98 0.82 
220 1.48 
270 2.70 


*Large temperature differences inside the reactor make calculation of 
actual velocities virtually impossible. For comparison purposes velocities 
are calculated assuming atmospheric pressure and 700°F. This figure 
includes both the dissociated ammonia for reduction and the séaling 
and cooling gases. The bed porosity averages 0.28. 


The Hydrofluorination Reactor 


The hydrofluorination reactor ft. high and 14-in. 
internal diameter. Its design capacity one mole UF, 
per hour (314 This reactor made entirely 
Inconel and equipped with three sets cooling coils 
using low-pressure steam the cooling medium. The 
steam discharged from these coils super-heated from 
the pick-up and used the reduction reactor for 
dhydration. The reactor also equipped with five 
blanket-gas spargers control hot spots the reactor. 
(See Figure 3). Hydrofluoric acid vapor 1100°F., fed 
into the bottom the reactor and flows countercurrently 
the pellets. The acid heated first steam jacketed 
carbon steel tank, which vaporizes the acid and heats 
140°F. six-pass Monel Brown Fintube raises the 
temperature the acid 275°F. and electric 
furnace with Inconel coil heats the acid 1100°F. The 


at 
1e 
SS ; 
re 
or 
1e 
e- 
age 
It, 
>SS : 


reactor equipped with five strategically placed vibrators 
maintain steady flow throughout the reactor. Five 
electric jackets with maximum heat output 
each are used for heating the reactor; they are controlled 
similarly those the reduction reactor. The temper- 
ature the bed the main reaction zone set 
1000-1100°F. Both reactors are controlled from central 
control panel. (See Figure 4). 

The stoichiometric excess acid from the reactor which 
runs about 10% passes through off-gas line along with 
the excess blanketing gas and scrubbed with aqueous 
caustic potash. The level the scrubber automatically 
controlled. order recover any entrained uranium 
carried over the off-gas line, the effluent filtered 
before being sent waste disposal. Alkalinity the 
scrubber liquor controlled continuous additions 
45% KOH small positive displacement pump. 

The rate product flow through each reactor con- 
trolled specially designed belt feeder located the 
bottom each reactor. The belt choke fed through 
seal leg. Table gives some typical data with regard 
the operation the reactor. 


TABLE 


OPERATING FLOWS AND PRESSURES THE 
REACTOR VARIOUS THROUGHPUTS 


in. 
Product Solids rate Gas velocity the reactor 
per hr. ft. per hr. ft. per sec.* off-gas line 
150 0.93 1.23 
190 1.17 1.86 
250 1.55 2.61 
290 1.80 3.35 


*Large temperature differences inside the reactor make calculation of 
actual velocities virtually impossible. For comparison velo- 
cities are calculated assuming atmospheric pressure and 700°F. This 
figure includes both HF vapor and sealing and cooling gases. It should 
be noted that velocities will fall after conversion of UO, to UF, 


because four moles of HF result in two moles of water being produced. 
The bed porosity averages 0.37. 


Conclusions 


Table illustrates the progress made the develop- 
ment acceptable nuclear grade green salt. 


TABLE 


PROGRESS PRODUCTION HIGH GRADE URANIUM 
TETRAFLUORIDE 


Length 

Start run run hours UF, 

Nov. 11, 1957 552 0.89 
Aug. 7,1958 696° 97.81 1.39 0.80 
0.50 0.75 


Mar. 1959 Continuous 98.75 


(a) Ammonium oxalate insoluble. 
(b) First run pilot plant. 
(c) Production plant product. 


Figure 4—Reactor control panel. The left hand side con- 
trols the reduction reactor and the right the 
tion reactor. 


may noted that development progressed the 
operation became much smoother, indicated the 
length runs made. 


From Table may observed that the production 
nuclear grade uranium tetrafluoride the moving bed 
process definitely feasible. offers many economic 
advantages over some other methods presently adopted— 
especially the use reagent chemicals, 
larly its adaptability for use smaller sized production 
units. 
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Evaluation Antoine Constants 
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The Antoine equation 


found adequate for representing and extrapolating 
vapor pressure data. method proposed for evaluating 
the constants which the approximate linear relationships 

the Clausius-Clapeyron plot (log vs. 1/T) 
employed. The proposed method provides much less 
involved procedure than the available calculation methods 
and does not involve the selection base-point 
which assumed free from error required 
the Thomson’s graphical method 

The constant may obtained evaluating the 
values log p/dt any two temperatures. From 
Equation (1) 


Therefore 


2.80 


220 


200 


280 300 320 340 3.60 
Figure 1—Vapor pressure benzene (Clausius-Clapeyron 
plot). 


and solve for However, not recommended 
determine the values the quantity log p/dt log 
vs. plot, slopes are difficult determine accurately 
such plot. Instead, advantage taken the fact that 
Clausius-Clapeyron plot, the resulting curve shows 
much less curvature and may considered linear 
over small temperature range. many instances, 
complete curve over moderate temperature range may 
approximately represented two straight lines. Two 
linear portions are taken from the plot and their slopes 
are used evaluate the constant 


slope log p/d (1/T) log p/dT) 
and 
0.5 0.5 
1 


The quantity taken the temperature corre- 
sponding the mid-point the linear section. 


2.80 
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Figure 2—Vapor pressure benzene (Antoine plot). 
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Therefore 
and 


220.0 


This method obtaining the constant believed 
simpler than the Thomson’s graphical method actual 
applications. especially useful when the vapor pressure 
data are only moderate precision. The value 
used construct Fig. where log plotted against 
1/(C 220.0), resulting good straight line which can 
conveniently extrapolated. The slope this straight 
line (—B). The constant may evaluated Equa- 
tion (1). The constants and may also solved 
means the equations 


and 


where the number points employed 
Figure The values and obtained are 


1204.6 


The vapor pressure data are used 
illustrate the procedure. Figure illustrates the plot 
log vs. 1/T, which six points are employed. Two 
linear sections are taken the plot corresponding the 
first two points and the last two points listed Table 


The calculated vapor pressure data are compared 
with the values taken from the Perry’s Handbook. 
evident that the agreement very satisfactory. 


TABLE 


COMPARISON CALCULATED VAPOR PRESSURE WITH 
LITERATURE VALUES 


6.8948 

7.6 1.6024 40.0 
15.4 1.7778 60.0 
26.1 100 2.0002 100.0 
42.2 200 2.3008 199.9 
60.6 400 2.6020 399.9 
80.1 760 2.8810 +0.3 

Nomenclature 


vapor pressure, mm. Hg. 
temperature, °C. 
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